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ABSTRACT 

We present numerical simulations of the axisymmetric accretion of a massive mag- 
netized plasma torus on a rotating black hole. We use a realistic equation of state, 
which takes into account neutrino cooling and energy loss due to nucleus dissociations. 
We simulated various magnetic field configurations and torus models, both optically 
thick and thin for neutrinos. It is shown that the neutrino cooling does not signifi- 
cantly change either the structure of the accretion flow or the total energy release of 
the system. The calculations evidence heating of the wind surrounding the collapsar 
by the shock waves generated at the jet- wind border. This mechanism can give rise to 
a hot corona around the binary system like SS433. 

Angular momentum of the accreting matter defines the time scale of the accretion. 
Due to the absence of the magnetic dynamo in our calculations, the initial strength 
and topology of the magnetic field determines magnetization of the black hole, jet 
formation properties and the total energy yield. We estimated the total energy trans- 
formed to jets as 1.3 x 10 52 ergs which was sufficient to explain hypernova explosions 
like GRB 980425 or GRB 030329. 

Key words: gamma-rays: bursts, methods: numerical, (magnetohydrodynamics) 
MHD, black hole physics. 



1 INTRODUCTION 

In spite of significant progress in recent years, the nature of 



gamma-ray bursts (GRB), discovered by |Klebesadel et al. 
( |1973| more than 30 years ago, is still enigmatic. Although 
the light curves and emission spectra of GRBs are very di- 
verse, they seem to split into two groups of possibly dif- 
ferent origin: long bursts (At > 2 s) with a softer spec- 
trum and short bursts (At < 2 s) with a harder spectrum 
Mazets fe Golenetskii||l981| |Kouveliotou et al.p993| |Fish- 



man &; Meegan 1995| ). The long GRBs are often believed 
to be associated with star-formation regions (Bloom et al. 



2002; |Fruchter et al.|2006]|Blinnikov et al.|2005| ). Only the"se 
regions can host massive stars that have astronomically very 
short lifetime and die soon after the birth. In fact, recent ob- 
servations have provided strong arguments in favour of the 
connection of GRBs with the deaths of massive stars. Light 
curves of many GRB optical afterglows show features inher- 
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ent in the supernovae events; moreover, several long GRBs 
have been firmly associated with particular supernovae, the 
most popular examples being GRB 980425 and SN 1998bw 



( Soffitta |1998) |Galama et aLl|1998| |Pian et al.||2000| ) . Even 
more convincing evidence exists in the case of the low red 
shift GRB 030329 (z = 0.1685; [Greiner et al.| p003)) and 



its associated supernova, SN 2003dh (iMatheson et al.|2003 



Hjorth et al.|2003|[Sokolov et al.|2003) . The spectra of these 
supernovae show exceptionally broad emission lines indicat- 
ing abnormally high velocity of the ejecta, typical of the 
recently proposed "hypernovae" class objects. 

The most popular model of the central engine of these 
sources is based on the "failed supernova" stellar collapse 
scenario implying that the iron core of the progenitor star 
forms a black hole ( |Woosley~|[T9 93 ) . If the progenitor does 
not rotate, its collapse is likely to happen 'silently' until 
all the star has been swallowed up by the black hole. If, 
however, the specific angular momentum of the equatorial 
part of the stellar envelope exceeds that of the last stable 
orbit of the black hole, then the collapse becomes highly 
anisotropic. While in the polar region it may proceed more 
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or less uninhibited, the equatorial layers form a dense and 
massive accretion disk. Then the collapse of the layers is de- 
layed, and the gravitational energy released in the disk can 
be very large and drive GRB outflows, predominantly in the 
polar directions where mass density of the accreting matter 
can be much lower (MacFadyen &; Woosley 1999). How- 



ever, the actual process responsible for the GRB outflows 
is not established and remains a subject of ongoing investi- 
gations. The main mechanisms proposed to explain GRB 
outflows are neutrino pair annihilation heating (Popham 



et al.|1999| |MacFadyen fe Woosley "[1 999 ; Alo y et al.|2000^ 



magnetized disk wind ( Blandford &; Payne||1982 Uzdensky 
&; MacFadyen ||2006 ), and magnetic braking of the central 



black hole rotation ( Blandford &; Znajek||1977| Barkov & 
Komissarov|2008a ) 



High-precision self- consistent models of disk dynamics 
and neutrino propagation are required in order to obtain 
reliable results in the neutrino-driven supernova explosion 
theory. By now only relatively crude calculations have been 
carried out, and they show that the neutrino heating need 
not have a dominant role. Birkl et al. (20071) studied the 



heating rate due to annihilation of neutrinos emitted by neu- 
trinospheres of various prescribed geometries and tempera- 
tures. The energy deposition rates obtained in this paper 
lie in the range (0.07 - 27) x 10 49 / erg/s, and the typical 
annihilation efficiency seems to be rather low, about 10 -3 . 
Neutrino heating from a geometrically thin standard accre- 
tion disk (Sha kura &; Sunyaevj 197*3 ) is calculated in a recent 
article by Zalamea & Beloborodov ( 2008 ) . It is shown that 



the process efficiency strongly depends on the rotation pa- 
rameters of the black hole and rapidly decreases with the 
distance as r -4-7 . Other aspects of the collapsar model were 



considered in Baushev & Chardonnet (2009). Shibata et al 



(2007) carried out general relativistic MHD simulations of 
accretion disks with the masses (0.1 — O.4)M . They found 
that the disk opacity for neutrinos was high, resulting in low 
(~ 0.01 — 0.02) efficiency of neutrino emission itself, as most 
neutrinos generated in the disk could not escape from it and 



accreted in the black hole. Nagataki et al. (2007) consid- 
ered both neutrino heating and cooling in their Newtonian 
simulations of collapsars. They concluded that the neutrino 
energy deposition was insufficient to drive GRB explosions 
and that the magnetic mechanism was more promising. 

In the last few years the role of magnetic field in driving 
the black hole accretion and relativistic outflows has been 
a subject of active investigations via numerical simulations, 
that produced numerous interesting and important results 
(|Proga et al.||2003| |De Villiers et"aLl|2003| [Hawley et al.| 
2006 ; |McKinney Gammie ||2004| |McKinney||2005[ 12006) 
Mc Kinney &; Narayan ||2007| [Beckwith et~. 2008a b;" [NaT] 



|gataki |2009||McKinney fe Blandford |2009| ). In these studies 
the initial distribution described a keplerian disk or equilib- 
rium torus threaded with a relatively weak poloidal field, 
whose lines followed the iso-density or iso-pressure contours 
of the disk. The disk accretion was found to be driven by 
magnetic stresses via development of magneto-rotational in- 
stability (MRI) (|Velikhov |1959| [Chandrasekhar |1960|[BaT- 



bus & Hawley 1991). In addition to the disk, the numerical 
solution considered two other generic structures - the mag- 
netized disk corona and the highly magnetized polar fun- 
nel that hosted the relativistic outflow from the black hole. 
These studies applied simple adiabatic equations of state 



and did not take into account effects of radiative cooling 
that may be small for some types of Active Galactic Nuclei, 
but not for disks of collapsing stars. 

Recently iShibata et al. (2007) have carried out two- 



dimensional general relativistic MHD simulations in Kerr 
metrics with a realistic equation of state and taking into 
account the neutrino cooling, but the physical time span of 
their computations was rather short, only ~ 0.06 s. 



In this article, we model the situation (similar to Beck- 



with et al. (2008a)) when a massive stellar core collapses and 
the outer envelop forms a massive accreting torus. In section 
2, we consider physical processes, in section 3 — initial con- 
ditions, in section 4 — results, section 5 is the discussion. 



2 PHYSICAL PROCESSES 
2.1 Equation of state. 

For the simulations we use the equation of state 

P g =P(p,T) = P (p)+n b kT+^, (1) 

where k is the Boltzmann constant, a is the radiation energy 
density constant, rib = p/m p is concentration of baryons, 
m p is the mass of proton, P is the pressure, p is the density 
and T is the temperature; Po(p) describes the ground level 
for the cold degenerate matter, the expression for Po(p) is 
an approximation of the table function from |Baym et al.| 
( |1971b|a| ); [Malone et al.] ( fl975l ) . 

The specific energy (per mass unit) was defined ther- 
modynamically as 



, n 3kT aT 4 
2m p p 



(2) 



The value of eo(p) was taken from Bisnovatyi-Kogan 



(2001J). This equation of state allows us to deal with high- 
density matter. A more accurate EOS is significantly more 
expensive in calculations, while the disregarding of the 
e + e~-pair component in the plasma, for instance, can give 
a relative error in definition of the temperature less than 0.3 
at high temperatures ( | Bisnovatyi-Kog an 2001). 



2.2 Neutrino cooling 

On the grounds discussed in the introduction we ignore ef- 
fects of neutrino heating in our simulations and take into 
account only the energy loss due to emission of neutrinos. 
Approximation formulas for the cooling rates due to e + e~ 
pair annihilation Q p , neutrino photo-production Q p h, and 
plasma mechanism Q p i are taken from|Schinder et ah ( 1987 ), 



the cooling rate due to URCA processes Q u — from 



et al. (1969), and the cooling rate due to synchrotron neu- 



trino emission Q s — from Bezchastnov et al. (1997) 



F u (p,T) = Q p + Q p h + Q P i + Qu + Q s - 



(3) 



The neutrino cooling is introduced by the source term in the 
energy-momentum equation 



-gSuf, 



(4) 



where u v is the plasma four- velocity, g is the determinant of 
the metric tensor, and S is the cooling rate as measured in 
the fluid frame. 



Accretion of a massive magnetized torus on a rotating black hole 3 





Figure 1. Initial conditions for the Low torus momentum (upper left panel), Quadrupole 1 (upper right panel), and Quadrupole 2 (lower 
panels) models. Initial distribution of the density logarithm and the initial magnetic field are shown by colours and lines, respectively. 



One of the aims of our research was the investigation of 
sensitivity of the torus accretion model to the neutrino cool- 
ing processes. We compare two cases: with no cooling at all 
and the optically thin case with neutrino cooling. The first 



instance is the most popular (Gammie et al.||2004| 


McKin- 


ney & Gammie |2004| |McKinney & Narayan |2007; 


Hawley 


et al. |2006| |Beckwith et al.| 2008b| and less realistic. The 



second instance well suits for a light torus with the mass 
of a few per cent of the solar. It could be realized when 
two neutron stars or a neutron star and a black hole merge. 
The intermediate optical depth case is considered in [Barkov 
(2008} and can be suitable for a collapsar model when the 



disk mass reaches several solar masses. 



3 SIMULATION SETUP AND THE 

DESCRIPTION OF CALCULATIONAL 
MODELS 

3.1 General consideration 

We make our calculation in the ideal relativistic MHD ap- 
proximation using an upwind conservative scheme that is 
based on a linear Riemann solver and the constrained trans- 
port method to evolve the magnetic field. The details of 
this numerical method and its testing results are expounded 
(1999); Komissarov (2004b); Komissarov 



Komissarov 



(2006) 



The gravitational attraction of the black hole is in- 
troduced via Kerr metric in the Kerr-Schild coordinates, 
{0, r, 0}. We set the black hole mass as Mbh — 10M© which 
corresponds to r g = 14.847 km (we define r g = GMbh / c 2 ). 
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Figure 2. Time evolution of the magnetic flux ^ through the sphere of radius r = 4.7r^ (left panel) and through the horizon (right 
panel) in the case of the Low torus momentum. The fluxes at the moments t = 4.96 • 10 -4 s., t = 0.0248 s., t = 0.0495 s., t = 0.0991 s., 
t = 0.346 s. are represented by solid, dashed, dot-dashed, dotted, and three-dots-dashed lines, respectively. 



The two-dimensional axisymmetric computational domain is 
(r in < r < r ou t,0 <6<tt), where r in = (1 + Vl - a 2 /2) r g , 
a = cJbh /GM bh is the dimensionless angular momentum 
of the black hole (Jbh is the black hole angular momen- 
tum), and Tout = 200r g ~ 2969 km. We adopt free-outflow 
as the outer and the inner boundary conditions. The inner 
boundary is located just inside the event horizon and adopts 
the free- flow boundary conditions. It is worthy of notice that 
the inner boundary is inside of the outer event horizon - this 
choice is possible since the horizon coordinate singularity is 
absent in the Kerr-Schild coordinates. The total mass within 
the domain is small (less than 25%) as compared with the 
mass of the black hole, which allows us to ignore its self- 
gravity. The grid is uniform over where it has 320 cells, 
and is almost uniform over log(r), where it has 459 cells, the 
linear size of each cell being the same in both directions. 
As the initial distribution we consider an equilibrium 



torus (jFishbone & Moncrief|[l976l jAbramowicz et al.|[l978 



Komissarov 



2006) which is a "torus" of plasma with the 



black hole in the centre. The value of the specific angular 
momentum lo determines the total effective potential which 
can be written as 




Figure 5. Time evolution of the maximum of the magnetic flux 
\I/ through the horizon for the models: Low torus momentum 
(dashed line), Neutrino cooling (solid line), and Quadrupole 2 
(dot-dashed-line) 



where 



and 



W(r, 0) 



■In 



: 9tct>gtct> — gug4>4> 



+ 2hg t( j ) + l gtt 



(5) 



(6) 



(7) 



Our numerical scheme cannot operate with vacuum. To 
avoid this difficulty, we have to introduce minimum pos- 



sible values of the density and pressure that limit how small 
they can be. We also check the magnetic field and do not 
allow it to drop below the value defined by the equations 

1 1 

pmin < -- — , and p m i n < — — -— . For the magnetization 

6 8ttc z 100 87r 

typical of the task considered these inequalities are usually 

satisfied with a margin of three orders. The magnetic field 

strength is limited by the pair production by the neutrino- 

antineutrino annihilation and the Compton up scattering 

pair production (Beski net al.|[l 992). 
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Figure 3. Distribution of sign(B r )(\B r |/10 16 G) 1//4 (B r is the radial component of the magnetic field). The upper panel represents the 
distribution at the moment t = 0.208 s. for the central part of the Low torus momentum model (left picture), general view of the Low 
torus momentum model (central picture), and Quadrupole 1 model (right picture). The lower panel represents the distribution for the 
Quadrupole 2 model at the moments t = 0.1288 s., t = 0.2278 s., and t = 0.4259 s. 




n t, & 



Figure 4. Time evolution of the magnetic flux ^ through the horizon for the Quadrupole 2 model. The left panel represents the time 
sections at t = 4.96 x 10 -4 s., t = 0.0248 s., t = 0.1238 s., and t = 0.4452 s. (solid, dashed, dotted, and three-dots-dished lines, 
respectively) 



3.2 Model choice 

The parameter we are most interested in is the energy release 
effectiveness 



Jo E t0t dt 
fo M BH c 2 dt 



(8) 



where E tot is the total energy flux at radius 180r^ and Mbh 
is the accretion rate on the horizon of the black hole. In 
order to test the effectiveness dependence on the accretion 
conditions, we introduce several models with various initial 
magnetic fields and other parameters. 

The initial magnetic field structure in the disk is mainly 
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determined by the dynamo effect. Dynamo in accretion disks 



was investigated by many authors (Torkelsson & Brander- 
burg||1994| |Arlt fe Riidiger||1999| |Bardou et al.||2001|) in 



the context of a-Q model. Depending on the dynamo pa- 
rameters, stationary dipole or quadrupole magnetic fields 
were obtained. Some models showed time oscillations. The 
calculations gave only a rough estimation of the magnetic 
field strength. To reproduce this simulations, we have cho- 
sen three different topologies of the initial magnetic field 
(similar to Beckwith et al. (2008a)). 

Due to a strong dependence of the toroidal field com- 
ponent on the poloidal one, we initially introduce a purely 
poloidal magnetic field. This field can be described by a vec- 
tor potential with a single nonzero component A^. In our 
calculations we use 



Ad, oc 



W(r, Of 
-W(r, 6>) 3 tanh( 

-^M) 3 tanh(^) 



e-Tv/2 

WQ 



dipole field 
Quadrupole 1 model 

Quadrupole 2 model 



(9) 

where we — 0.05 and w r = 0.1 are the parameters of the 
vector potential. We adopt a dipole initial magnetic field 
for Low torus momentum, High torus momentum, Neutrino 
cooling, Schwarzs child, and Low magnetized models. Model 
Quadrupole 1 possesses a quadrupole-like field generated by 
two dipoles: one is above and the second is below the equa- 
torial plane. Quadrupole 2 model also has a quadrupole field, 
but the two dipoles are situated closer and farther from the 
black hole, which reproduces time oscillating magnetic field 
in the accretion disk. The topology of magnetic field is pre- 
sented in Fig.^ In all the cases the initial field is normalized 
so that the maximum of the magnetic to gas pressure ratio 
(3 = Pm/Pg does not exceed 3-10 -2 , except for the Low mag- 
netized model, where f3 < 3 • 10 -4 . We have chosen the small 
value of /3 in order to avoid strong influence of the magnetic 
pressure on the initial hydrodynamical equilibrium of the 
configuration. 

The dimensionless angular momentum of the black hole 
is a = 0.9 for all the models except for the Schwarzschild 
one, where a = 0. The collapsing torus specific angular 
momentum and mass are equal to / = 2.8 • r g c = 1.247 • 
10 17 cm 2 sec -1 , Mtor — 2.55Mq for Low torus momentum, 
Neutrino cooling, Quadrupole 1, Quadrupole 2 models and 
I = 4.0-r^c = 1.782- 10 17 cm 2 sec' 1 , M tor = 1.87M© for High 
torus momentum, Schwarzschild, Low magnetized models. In 
all the models, except for Neutrino cooling, we neglect the 
neutrino cooling. The summary of the model parameters see 
in tab. 1. 



stability (MRI) defines the accretion rate and could be the 
source of the magnetic dynamo mechanism in the disk. Due 
to the 2D approach MRI cannot lead to the magnetic dy- 
namo appearance in our calculations. It is easy to see that 
by the magnetic flux 



Jo 



B p dS 



(10) 



decreases with time (Fig. [5] Low torus momentum model). 
First of all, we investigate the influence on our system of 
various magnetic field topologies. 

A dipole-type Low torus momentum model exhibits a 
fast MRI growth in the disk. At first, we can see a powerful 
burst of energy extraction. It seems to be an effect of switch- 
ing, after a while the flux becomes stationary. The accretion 
runs actively. Magnetosphere of the black hole is dipole-like 
(see Fig. [5]). The magnetic flux through the horizon is of the 
same order as the initial one; the effectiveness of the total 
energy extraction 77 is ~ 7.4 x 10 -4 . 

Quadrupole 1 model shows a slower growth of MRI. The 
magnetosphere of the black hole is monopole-like (see Fig. [3]) 
and does not significantly change with time. The magnetic 
flux through the horizon is ten times reduced, compare to 
the initial one that results in a low energy extraction rate: 
f] is 1.1 x 10" 5 . 

The magnetic field structure of the Quadrupole 2 model 
actually represents a configuration of two current frames on 
the different radii in the equatorial plane. In this case, the 
magnetosphere is dipole-like, but time variant. One can see 
on Fig.[3]that initially the magnetic lines near the black hole 
are predominantly oriented upward. As the inner currant 
frame is falling under the horizon, the magnetic field lines 
get the opposite direction. The mean magnetic flux is also 
ten times less than the initial one. The efficiency coefficient 
j] of the jet is 2.2- 10~ 5 . 

In our calculations the accretion rate to the black hole 
is proportional to the maximum value of the magnetic flux 
module in the torus. On the other hand, the energy release 
in both the quadrupole cases drops down dramatically (up 
to 250 times, see tab. 2). 

In the case of the Low magnetized model, the MRI grows 
very slowly. Accretion starts after a long (0.6 sec) phase of 
the field linear amplification in the torus up to (P m /P g ) ~ 1. 
When the accretion starts, the system takes on the standard 
appearance: a jet, an accretion disk and the wind from it. 
The energy release and the accretion rate are very low (see 
tab. 2). The efficiency r\ of the wind is2.0xl0 -5 and of the 
jet is 4.9 x 10" 7 . 



4 RESULTS 

All the models (except the Schwarzschild) after the first 
relaxation gave rise to the same standard configuration: a 
highly magnetized jet, a thick disk and the wind from it. 
Let us discuss the influence of various parameters of the 
system on the accretion rate and energy yield. 

4.1 Magnetic field topology 

Poloidal magnetic field plays a crucial role in the evolution 
and dynamics of the accretion disks. Magneto-rotational in- 



4.2 Angular momenta of the torus and of the 
black hole 

In the models Low torus momentum, Quadrupole 1, and 
Quadrupole 2 we put the torus angular momentum / = 
2.8 • r g c which is very close to the limit configuration, and 
even small perturbations lead to accretion. If the initial mo- 
mentum is augmented to / = 4.0 -r g c (High torus momentum 
model), it decreases the accretion rate thirty times but the 
energy release does not drop so strongly (see tab. 2), since a 
significant part of the field accretes on the black hole increas- 
ing the magnetic flux near it. As a result, the effectiveness 
of energy extraction in this case is high 77 = 3.0 x 10 -3 . 
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Table 1. The list of model initial parameters. 



model name 


l/(r g c) 


a 


Mtor 


Neutrino cooling 


Magnetic field type 


max(Pm/P g ) 


Quadrupole 1 


2.8 


0.9 


2.55 


No 


Up-down Quadrupole 


0.03 


Quadrupole 2 


2.8 


0.9 


2.55 


No 


Farther-closer Quadrupole 


0.03 


Low torus momentum 


2.8 


0.9 


2.55 


No 


Dipole 


0.03 


Neutrino cooling 


2.8 


0.9 


2.55 


Yes 


Dipole 


0.03 


High torus momentum 


4.0 


0.9 


1.87 


No 


Dipole 


0.03 


Schwarzs child 


4.0 


0.0 


1.87 


No 


Dipole 


0.03 


Low magnetized 


4.0 


0.9 


1.87 


No 


Dipole 


0.0003 




Figure 6. Dependence on 9 and time of the matter (E MA , left panels) and electromagnetic (E EM , right panels) energy fluxes per unit 
angle at the radius R = 180r g for the Schwarzschild (upper row) and High torus momentum (lower row) models. 



Another important parameter is the angular momen- 
tum of the black hole a. Detailed investigations of its influ- 
ence on the jet formation has been mad e m[G ammie et al.| 



(2004); McKinney fc Gammie (2004); McKinney (2005); 



Hawley et al. (2006). The structure of the solution signifi- 



cantly changes if a — (Schwarzschild model, see tab. 2). 
The central highly magnetized jet does not appear at all, 
the energy is released through the subrelativistic magne- 
tized wind. The total energy flux of the system jet- wind is 
almost the same for the Schwarzschild and the High torus 
momentum models, whereas the accretion rate in the former 
case is 27 times higher. 

Let us consider the time-space evolution of the energy 



flux at the distance R = 180 • r g . The total flux can be di- 
vided in two parts: electromagnetic E EM and carried out by 
matter E MA . In the case of High torus momentum (Fig.pl 
the lower row), it is easy to see that the flux is predom- 
inantly electromagnetic in the jets and matter-dominated 
in the wind. The maximal electromagnetic flux propagates 
near (not on) the axes. The maxima of the matter energy 
flux are strongly correlated with the electromagnetic max- 
ima, as the main part of matter energy flux is driven by 
shock waves appearing at the boundary between the rela- 
tivistic jets and the non-relativistic wind (the lines going 
up from the jet to the equatorial region). The jets very ef- 
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fectively heat the surrounding wind and can produce a hot 
corona around the disk. 

In the case of Schwarzs child model, the situation is dif- 
ferent (fig(6j the upper row): electromagnetic and matter en- 
ergy fluxes have much wider angular distribution. A highly 
magnetized zone does not appear, and we can see a mat- 
ter dominated magnetized outflow. Magnetic braking of the 
black hole produces the wind that carries out the magnetic 
field. 



4.3 Neutrino cooling 

In the initial torus configuration radiation and degenerated 
matter produce 80% and 20% of the total pressure, respec- 
tively. Therefore, neutrino cooling becomes very important. 

Intensive neutrino cooling in the Neutrino cooling model 
(where the system is not opaque for neutrino propagation) 
leads to a significant disbalance in the initial distribution. 
The originally stable torus collapses to a new configuration 
with negligible thermal pressure. One can see the influence 
of neutrino cooling on the structure of the torus in Fig. [7] 
The torus without cooling keeps the initial size and shape. 
The cooling leads to a collapse of the initial configuration. 
The maximum density increases up to ten times and reaches 
the value 10 13 g/cm 3 . 

The structure of the accretion flow does not change sig- 
nificantly: the Neutrino cooling model flow is similar to High 
torus momentum (Fig. [7| and thinner than in Low torus 
momentum model. The Low torus momentum model has a 
much smaller and less stable highly magnetized jet region 
than in the Neutrino cooling case (Fig. [8]). Despite the fact 
that the Neutrino cooling accretion rate is two times lower 
than the Low torus momentum one (tab. 2), the energy re- 
lease and the magnetic flux through the horizon are even 
slightly higher in the Neutrino cooling model (Fig. [5}. So 
Neutrino cooling and High torus momentum models have 
the same effectiveness coefficient 77 ~ 0.003 that is 2 — 4 
times greater than for the Low torus momentum model. 

At the inception of the calculations the central region 
of the torus rapidly cools down and falls to the potential 
hollow, while the outer parts do not cool at first. Inten- 
sive neutrino cooling leads to a disbalance in the accretion 
torus. Due to asymmetry of the potential well, uncompen- 
sated gravitational forces give rise to oscillations (Fig. |9j 
similar results were obtained by Zanotti et al. (2005); Mon-| 
tero et aLj ( |2007| ). So the cooling produces intensive ra- 
diative shock waves, that (as well as the poloidal magnetic 
field) strongly increases the effective viscosity. When we ne- 
glect the cooling, the oscillations do not appear, the mean 
radius increases due to the strong wind outflow. In the op- 
tically thin case, the oscillations appear with the amplitude 
Ar m /r m 0.05, where r m f M \r\dm/ f M dm is the mean 
mass radius, and quickly degrade due to the high viscosity. 
We can estimate the relaxation time as t re iax = t osc x Q, 
where Q = |(Arf + ArJ^^Arf - Ar£i)| ~ 3.8 and 
Ar™ is the maximal deviation from the equilibrium of the 
i-th oscillation. In our case t re iax ~ 0.027 sec. The same os- 
cillations were obtained in the case of intermediate neutrino 
optical depth ( |Barkov ||2008| ). 




0.04 0,06 
time, sec 

Figure 9. Mean radius r m evolution for Low torus momentum 
(dashed line) and Neutrino cooling (solid line) models. 



5 DISCUSSION 

Topology and strength of the magnetic field play a crucial 
role in the process of accretion on a black hole. In our simu- 
lations the energy release effectiveness varies over two orders 
of magnitude when we change the magnetic field topology 
from the dipole (Neutrino cooling, High torus momentum 
models) to the quadrupole one( Quadrupole i, Quadrupole 2 
models). Due to the 2D approach we cannot consider the dy- 
namo effect. However, we can see temporal increases of the 
magnetic field intensity owing to the MRI-driven turbulence 
growth and magnetic line twist. 

The Blandford-Znajek mechanism plays a leading role 
in accreting black hole energy release. For a force- free 
monopole magnetosphere the Blandford-Znajek power is 
given by 



Ebz 



1^ { Q^f 
6c V 47 



(11) 



We assume Q = Qh/2, were Qh = a/2(l + y/l — a 2 ) is the 
angular velocity of the black hole and ^ is the magnetic flux 
( 10 ) threading the black hole. 

Our estimations of the energy release effectiveness differ 



by a factor of 3 from McKinney (2005) and by a factor of 
up to 25 from Barkov & Komissarov ( 2008b ) ; | Komissarov| 



&; Barkov ( 2009 ) . We explain it by distinctions between the 
initial magnetic fields: in Komissarov & Barkov ( 2009| , for 
instance, the magnetic field flux was 4.8 times higher. 

Neutrino cooling processes significantly suppress the 
wind from the accretion disk, but do not diminish apprecia- 
bly the energy release; moreover, sometimes they can sta- 
bilize the outflow. The neutrino driven oscillations (Fig. [9]) 
look like an artificial effect of turning on the initial condi- 
tions: more realistic conditions like MacFadyen et al. ( 1999); 
Barkov & Komissarov (2008al do not lead to significant ra- 



dial oscillations. But even if the effect is real, the oscillations 
( |MacFadyen et al.|1999) can scarcely be observed by gravi- 
tational wave detectors. 

Highly magnetized jets appear in all the models con- 
sidered, where the black hole possesses angular momentum. 
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Figure 7. Distribution of log 10 (p) and magnetic lines. Left panel represents the initial distribution. Central and right panels represent 
the distributions after 0.2075 s. of evolution for the Low torus momentum and Neutrino cooling models are shown on the central and 
right panel, respectively. 





Figure 8. Lorenz factor distribution for t = 0.2075 s. Left panel represents the Low torus momentum model with no cooling, the right 
one — Neutrino cooling model. 



Table 2. Main results. 



model 


I 


a 


M BH 


M w 




-hMA 
^51 




V 


Low torus momentum 


2.8 


0.9 


1.3937 


0.0189 


0.777 


1.076 


1.8530 


7.4 x 10~ 4 


Neutrino cooling 


2.8 


0.9 


0.8588 


0.0134 


1.5024 


1.4014 


2.9039 


1.9 x 10- 3 


Quadrupole 1 


2.8 


0.9 


0.3976 


0.00023 


0.0022 


0.0058 


0.0080 


1.1 x 10- 5 


Quadrupole 2 


2.8 


0.9 


0.2126 


0.00034 


0.0019 


0.0065 


0.0083 


2.2 x 10- 5 


High torus momentum 


4.0 


0.9 


0.053 


0.0026 


0.1283 


0.1541 


0.2824 


3.0 x lO- 3 


Schwarzschild model 


4.0 


0.0 


1.411 


0.0086 


0.0216 


0.2320 


0.2536 


1.0 x lO" 4 


Low magnetized 


4.0 


0.9 


0.0023 


5.5 x 10- 6 


2.0 x 10~ 6 


8.16 x 10~ 5 


8.36 x 10- 5 


2.0/0.049 x 10 



Here I is the dimensionless angular momentum of the accreting matter; Mbh — accretion rate Mqs -1 ; M w - wind mass loss rate 
Mq s _1 at radius 180r^; E^[ M , E^f A , are electromagnetic energy, matter energy, and the total energy fluxes at radius 180r^, 

10 51 erg s _1 ; n — accretion effectiveness ([SJ (for the Low magnetized model we present separately the effectiveness of the 
electromagnetic energy extraction). Notice that in the case of the Low magnetized model the accretion starts only after 0.75 s., while 

the calculation ends at 0.99 s. 
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The obtained parameters of the plasma acceleration in the 
jets (Fig . |8| are in good agreem ent with 
(2007a|); iBarkov & Komissarovl d2008bb; 



(2009). 
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The boundary between the jet and the wind is unsta- 
ble (the perturbation growth is most pronounced in Fig. [8] 
Neutrino cooling case, in the south hemisphere in the ra- 
dius interval 100 — 150r g . [McKinney] ( |2006| reported about 
a similar effect). The jet speed is much higher than the sound 
speed in the wind, and the boundary perturbations generate 
strong shock waves propagating to the wind region (Figj6|. 
The shocks effectively heat the wind and can produce a hot 
corona around the accretion disks like in the SS433 system 
( |Cherepashchuk et al.|2006 ). 

If we adopt the effectiveness coefficient 77 ~ 0.003 of 
the models Neutrino cooling and High torus momentum, it 
allows us to estimate the total energy release of the torus 
accretion as E to t ~ r/MtorusC 2 w 1.3 • 10 52 erg. This value 
of 77 however, need not to be universal: as we saw, other 
magnetic field topologies can have a significantly different 
efficiency. 

The scenario considered in this article is a possible al- 



ternative to the magnetar-driven hypernova scheme ( Komis- 



sarov &; Barkov |2007b| ). Both the scenarios predict the ap- 
pearance of collimated jets energetic enough to explain GRB 
as well as hypernovae phenomena. Accurate calculations are 
needed to find characteristic features and distinguish these 
two schemes. By now we can say that the chemical com- 
position of matter in the jets should be different. In the 
collapsar scenario we expect a highly Poynting-dominated 
jet (efieid/ematt ~ 1000, where ef ie id and Cmatt are the field 
and the matter energy flux densities, respectively), mainly 
made up of e + e~ -pairs, while in the magnetar-driven sce- 
nario the jet is baryon-dominated and intermediately mag- 
netized (€field/ Cmatt ~ 1 - 100). 

The effect of the magnetic field sign change found in 
the jet of the Quadrupole 2 model, may be the key to un- 
derstanding of the "surprising evolution of the parsec-scale 
Faraday rotation gradients in the jet of the BL Lac object 
B1802+784" ( |Mahmud et al. ||2009| ). Complicated topology 
of the magnetic field in the accretion disk can lead to varia- 
tions of the magnetic field direction in the jet. In this case, 
the jet constitutes a chain of regions with various magnetic 
field orientations. 



6 CONCLUSION 

The accretion rate and energy release are governed by four 
main parameters: magnetic field structure, angular momen- 
tum of the black hole, specific angular momentum of the 
accreting matter, accretion disk mass or accreting part of 
the star envelop. To have a long-time accretion, angular mo- 
mentum of the accreting matter should be high enough to 
form an accretion disk around the black hole. As it has been 
demonstrated by the Schwarzs child model, if the black hole 
has deficient angular momentum, the powerful magnetized 
jet does not appear at all, which totally changes the picture 
of the accretion. Moreover, a comparison of the Low and 
High torus momentum models shows us that a respectively 
small variation of the angular momentum of the accreting 
matter leads to a change of the accretion efficiency. Neutrino 



losses do not prevent energy release and can even stabilize 
the magnetized jet outflow. 

The least understood and a very important factor defin- 
ing the accretion efficiency is the disk magnetic field struc- 
ture. As we could see from the calculations, a mere replacing 
of the dipole field by the quadrupole one decreases the en- 
ergy release by two orders of magnitude. The magnetic flux 



^ governs the energy extraction from ergosphere (11), while 



the black hole magnetosphere is formed and supported by 
the disk. In our calculations we preset the initial magnetic 
field in the torus. It is clear, however, that self-consistent 
numerical calculations of the magnetic dynamo in the disk 
are necessary in order to elucidate the problem. 

Nevertheless, our calculations show that a very inten- 
sive electromagnetic and matter energy flux can be evolved 
during massive torus accretion on a rotating black hole 
with mass ~ IOMq. The total energy budget is quite suf- 
ficient to explain hypernova explosions like GRB 980425 or 
GRB 030329. The magnetized matter accretion on a rotat- 
ing black hole is a very promising model of the central engine 
of these phenomena. 

In addition, we have found that the instability of the 
boundary between the jet and the wind generates shock 
waves propagating to the wind region and heating it. This 
effect can give rise to a hot corona around the binary system 
similar to SS433. 
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